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The electronic structure of six low-lying electronic states of scandium hydride, X a®A, 

fe^n, A^A c 3E+, and B^n, is studied using multi-reference configuration interaction as a 
function of bond length. Diagonal and off-diagonal dipole moment, spin-orbit coupling and 
electronic angular momentum curves are also computed. The results are benchmarked against 
experimental measurements and calculations on atomic scandium. The resulting curves are 
used to compute a line list of molecular ro-vibronic transitions for "^^ScH. 
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1. Introduction 


Scandium hydride was first identified experimentally by Smith [l|, who recorded 
the absorption spectra of various transition metal hydrides (ScH, TiH, VH, NiH, 
CoH and deuterated isotopologues) in the region 17700 to 18300 cm“^. No detailed 
analysis of the spectrum was reported but it was remarked that a triplet ground 
state was expected. Later studies showed that the ground state of ScH is actually 
S+j with a low-lying '^A. Early studies using restricted open-shell Hartree-Fock 
(ROHE) [3, S, generalised valence bond theory [3] and empirically-fitted pseu¬ 
dopotentials a] all incorrectly predicted a ground state, with the S"*" generally 
lying about 2000 cm“^ higher up. These studies considered the six electronic terms 
correlating upon dissociation with ground state atoms (dissociation channel la¬ 
belled 1 in table m leading to ^’^H and ^’^A); it was remarked 0] that the 
bonding of the molecular terms other than is due to the Sc(4s) and H(ls) 
orbitals, while the scandium 3d orbitals are relatively unaffected with respect to 
the atomic state [3]- On the other hand the bonding of was put down to spd 

bonding [3, 0]- The different bonding character of the term is probably one 

of the reasons why an extensive treatment of electron correlation is necessary to 
obtain the right ordering of the electronic terms. Bauschlicher and Walch [3] were 
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the first to correctly predict S"*" lying below A by performing full-valence multi¬ 
configuration self-consistent-field (MCSCF) calculations. Jeung and Koutecky Q] 
studied the same six electronic terms using pseudopotentials and truncated MRCI 
and confirmed a ground term close to equilibrium (re ~ 3.4 ag), although 
for longer bond lengths becomes lower in energy. Note that all these studies 
kept the scandium outer-core 3s3p electrons uncorrelated and often did not include 
relativistic corrections. 


Anglada et al produced a series of papers [lfll4l3l| studying in great detail ScH and 
ScH'*' using multi-reference configuration interaction (MRCI); in particular their 
hnal paper (131] constitutes the most complete theoretical study of ScH currently 
available. These authors confirmed that is the ground state when correlation 
effects are included; they also found that correlation of the Sc{3s3p) semi-core 
electrons leads to large energy shifts, strongly stabilising the term with respect 
to the others and swapping the order of some of the excited states. They used basis 
sets similar in size to cc-pVTZ. 

More recent theoretical studies considered ScH in the context of calibration stud¬ 
ies of transition metal molecules using density function theory [l^ . 15|. but they 
focused on equilibrium properties of very few terms and are of little relevance for 
us. An exception is the very recent study by Hubert et al 16( in which a detailed 
study of the ground and of two excited terms ^’^A around equilibrium was 
presented and a modification of the coupled cluster method called general active 
space coupled cluster (GAS-CC) was used. 

The only theoretical dipole moment data available for ScH are those by Anglada 


131] and by Chong et al 17]. 


et al 

Experimentally a study by Bernard et al 1^ reported three new bands ascribed 
to ScH and ScD in the region 11600 to 12700 cm“^, but no detailed analysis or 
assignments were made due to the limited resolution and complexity of the spectra. 
More recently Ram and Bernath [l^ . 20] reported two detailed emission spectra 
analyses for ScH and ScD. In these studies they reported on singlet-singlet bands 
in regions from 5400 to 20500 cm“^ assigned to transitions between 8 electronic 
terms, namely A^A, 0^2+, D^H, E^A, and G^H. Two 

additional strong bands near 11620 and 12290 cm“^ and two weaker bands near 
12660 and 16845 cm“^ were recorded but only incompletely analysed; the band near 
11620 cm“^ was conjectured to be due to a transition to the low-lying ^A term 


211 ] reported more recently a detailed experimental 


from a term. Le and Steimle 
study of the X ^H band around 16850 cm“^, where also the electric dipole 

moments of ScH in its and states were obtained using optical Stark 

spectroscopy. Very recently, Mukund et al [2^ reported the observation of ScH 
emission bands at about 17900 cm“^, ascribed to the g - a ^A triplet-triplet 
transitions. 

This paper focuses on the six low-energy electronic states dissociating to ground 
state Sc and H atoms. Of the various experimentally observed bands 1914211] only 
X - B is considered in the present study, although experiment is used for the 
empirical refinement of the potential energy curves of the singlet terms X A ^ A 
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and as well. 

As part of the ExoMol project 23|], whose aim is to produce comprehensive 
line lists for hot, astrophysically-important molecules, we have been constructing 
rovibrational and rovibronic line lists for a number of diatomic species 24l-l28l|. 
However, so far none of these have contained a transition metal (TM). The rich¬ 
ness of the spectrum of TM-contain ing diatomics makes their opacity particularly 
important for astrophysical studies 29| but treating their rovibronic spectrum ab 
initio is very challenging. Scandium hydride is the lightest TM molecule and for 
this reason constitutes a useful benchmark system for theoretical studies of such 
systems. 

Scandium hydride has received comparatively little attention with respect to 
other transition metal hydrides such as FeH or NiH, in all probability because of 
the low abundance of scandium. Scandium is in fact the rarest of fourth-period 
transition metals (Sc-Zn) in the solar system 30(, although it is more abundant 
than all heavier elements starting from the fifth period. The study presented here 
on ScH constitutes a first step in the ab initio calculation of ro-vibronic spectra of 
TM-containing molecule. We perform a series of ab initio calculations on both the 
scandium atom and ScH and use these to produce a line list of scandium hydride 
line positions and intensities which is reasonably complete in the region up to 
12 000 cm“^. 


2. Atomic Scandium 


As a preliminary test we studied in some detail the scandium atom using com¬ 
plete active space self consistent field (CASSCF) and internally-contracted MRCI 
calculations using the program Molpro 311]. 

We collected in tabled] reference energy levels of the scandium atom up to about 
20 000 cm“^, along with the ScH molecular terms correlating adiabatically with the 
various atomic states; this information serves as an indication of which molecular 
terms are expected to be low-lying; the rationale is that for transition-metal hy¬ 
drides the hydrogen atom constitutes a relatively small perturbation of the atomic 
energy levels, so that molecular terms correlating with high-energy atomic prod¬ 
ucts should be high-lying too. The lowest dissociation channel is separated by about 
11 500 cm“^ from others, and in fact the six electronic terms correlating to it are 
lowest-lying and most theoretical studies concentrated on them. The dissociation 
channels 1 to 10 reported in table[T]lead altogether to 60 electronic molecular terms; 
that is, there are 60 energy curves in absence of spin-orbit splitting, which become 
155 in presence of spin-orbit splitting. Of these, only eight have been characterised 
experimentally 20j . Very probably many of the remaining molecular energy curves 
are repulsive (i.e., have no minimum) and therefore do not concern us here. In 
any case great complexity and strong perturbations in the observed spectra are 
expected. This is typical of open-shell transition metal molecules. 

In the following we consider the excitation energies from the ground term 
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Table 1. Energy levels for scandium atom up to 20 200 cm~^. The {E) are term-averaged energies, com puted 
by {E) — ’ experimental energies for the levels Ej where taken from the NIST website I 32 I 1 . n/ — 

min(2S -I- 1,2L -|- 1) is the number of fine-structure components of a given terms due to spin-orbit interaction; 

— {2S -|- 1)(2L -|- 1) is the total degeneracy of the term (number of microstates). The ^’s are effective spin- 
orbit coupling constants, such that the spin-orbit splittings for each term are best reproduced by the expression 
Ej — Eq -hi) — L{L -h 1) — S{S -h l)]/2, J = \L — S'!, ■ ■ ■ , L -\- S. The last column lists the molecular 

terms for ScH correlating at dissociation with the given Sc atomic term plus a ground state ^S hydrogen atom. 


# 

config. 

term 

(E) / cm ^ 

Uf 

nfj, 

^ / cm ^ 

molecular terms 

1 

SdHs^ 

2d 

0.0 

2 

10 

67.3 

h3[s+,n,A] 

2 

3d^4s^ 

4f 

11 509.1 

4 

28 

15.0 

A,$] 

3 

3d^4s^ 

2F 

14 891.3 

2 

14 

33.1 

i’3[s-,n,A,ci>] 

4 

Sd^As^Ap^ 

4 po 

15 775.8 

4 

28 

33.6 

3’'5[s+,n, A,$] 

5 

?,d}As^Ap^ 

4do 

16 031.0 

4 

20 

25.2 

3.5[s-,n,A] 

6 

?,d}As^Ap^ 

2D° 

15 951.4 

2 

10 

-29.7 

h3[s-,n,A] 

7 

2,d?As^ 

2d 

16 916.7 

2 

10 

-5.0 

h3[s+,n,A] 

8 

3^245! 

4 p 

17 175.2 

3 

12 

20.1 

3.5[s-,n] 

9 

Sd^As^Ap^ 

4 po 

18 440.6 

3 

12 

15.0 

3’'5[s+,n] 

10 

4s‘^4p^ 

2 po 

18 706.5 

2 

6 

96.5 

i’3[s+,n] 


to the two lowest-energy excited terms, namely and ^F. We expect that errors 
of computed molecular excitation energies are comparable with the corresponding 
error in the atomic case. 

The most accurate study of (neutral or singly-ionized) transition metal atoms 
including scandium is due to Balabanov and Peterson [33, ISJ] . These authors used 
coupled cluster up to CCSDTQ, relativistic corrections based on the Douglas-Kroll- 
Hess (DKH) hamiltonian, included core correlation and developed the largest basis 
sets available for transition metals. For scandium —)■ '^F excitation energy the 
best theoretical coupled cluster result is higher than the experimental one by about 
115 cm“^. The corresponding result using ACPF (a multi reference method very 
close to MRCI) in conjunction with the full-valence reference space is too low by 
about 110 cm“^; using a larger reference space including a further set of diffuse 
d functions (which are thought to be necessary for describing the late transition 
metals Fe-Cu) leads to a worse agreement with experiment of about 190 cm“^. 
Other recent studies of transition metal atoms excitation energies including scan¬ 
dium were performed by Raab and Roos (35(| and Mayhall et al [36(. Raab and 
Roos 351 computed the ^F excitation energy with CCSD(T) and CASPT2 


using the DKH hamiltonian for relativistic effects and the ANO-RCC basis set 
(similar in size to aug-cc-pCVQZ). Both CCSD(T) and CASPT2 frozen core val¬ 
ues agree with experiment to about 250 cm“^, but allowing for core correlation 
worsens somewhat the agreement to about 500 cm“^. Mayhall et al 36|] also used 
core-correlated CCSD(T) with the G3Large basis set (similar in size to aug-cc- 
pCVTZ) and reported an agreement of 250 cm“^ without inclusion of relativistic 
effects and of about 1200 cm“^ when relativistic effected were included. 

Table [T] gives some indicative result for both the ^F and —>■ ^F transi¬ 

tions performed in this study using MRCI and the full valence reference space. Our 
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Table 2. Electronic term excitation energies for scandium atom (this work). All calculations used the full-valence 
(3-electron, 9-orbital) complete active space comprising the 3d4s4p orbitals. Orbitals are state-averaged over the 
three electronic terms considered. MRCI-I-Q are Davidson-corrected energies (relaxed reference). Calculated values 
are reported as (reference — calculated), and reference energies are taken from the column labelled ‘(E)’ of table 
m All quantities are in cm~^. 


Method 

reference 

energies= 

^ ^F 2d ^ 2p 

11 509.1 14 891.3 

ref - calc 

CASSCF 

3z 

28.6 

64.0 

CASSCF 

4z 

49.2 

86.0 

MRCI/frz core 

3z 

-1187.9 

-421.8 

MRCI/frz core 

4z 

-1084.6 

-265.3 

MRCI/core corr 

wc3z 

573.3 

541.1 

MRCI/core corr 

wc4z 

1072.1 

1056.0 

MRCI-l-Q/core corr 

wc4z 

2080.0 

2039.8 

MRCI-|-Q/core corr/DKH4 

wc4z-DK 

813.6 

707.1 


best results for both transitions are too small on the average by about 750 cm“^ 
with respect to the results by Balabanov and Peterson; our errors are larger prob¬ 
ably because we performed a state-averaged calculation at the CASSCF level and 
also because of the smaller basis sets used. A striking consideration is that the 
non-relativistic CASSCF excitation energies are extremely good, a fact which can 
only be due to fortuitous cancellation of errors. Overall our results and the analysis 
of the literature show that it is difficult to get excitation energies correct to better 
than about 500 cm“^, and that good agreement with experiment can be often due to 
cancellation effects. We also observed that the relativistic contribution to excitation 
energies shows relatively large variations of the order of 500 cm“^ depending on the 
levels of electron correlation (CASSCF, MRCI valence only or core-correlated) and 
on using the mass-velocity one-electron Darwin (MVDl) rather than the Douglas- 
Kroll-Hess (DKH) Hamiltonian. 

We also computed atomic spin-orbit splitting constants using CASSCF and 
MRCI wave functions as implemented in MOLPRO. Results are collected in ta¬ 
ble [3l Spin orbit splitting constants show weak sensitivity to the size of the basis 
set and already with the smallest 2z basis set are converged within 1 cm“^. The 
dependence on the electron correlation treatment is also weak, with the ground 
term being the most sensitive. Going from CASSCF to frozen-core MRCI in¬ 
creases ^(^D) by -|-18 cm“^ but reduces ^(^F) and ^(^F) by only 1.5 and 0.9 cm“^ 
respectively. Correlating the (3s3p) outer core increases by 5 cm“^, ^(^F) by 

1.8 cm“^ and ^(^F) by 2.6 cm“^. With respect to the experimentally-derived values 
we do not observe a clear pattern of convergence with respect to the level of theory 
used, and the simplest CASSCF/2z values agree with experiment practically as 
well as the core-correlated, large basis set MRCI ones. 
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Table 3. Calculated spin-orbit constants ^ for scandium atom. The column labelled ‘obs.’ are experimentally 
derived values from tabl elTl All values are in cm~^. _ 

CASSCF 


transition 

Obs. 

2z 

3z 

4z 

5z 

C(2D) 

67.3 

57.5 

57.6 

57.9 

57.9 


15.0 

18.5 

18.5 

18.6 

18.6 

|(2f) 

33.1 

37.2 

37.2 

37.4 

37.4 



MRCI (frz. core) 


transition 


2z 

3z 

4z 

5z 


67.3 

74.7 

75.6 

76.0 

76.0 

CCf) 

15.0 

16.8 

17.0 

17.1 

17.1 


33.1 

35.9 

36.4 

36.5 

36.6 



MRCI (core 

correlated) 

transition 



wc3z 

wc4z 

wc5z 

C(2d) 

67.3 


80.3 

81.2 

81.3 

l('F) 

15.0 


18.7 

18.9 

19.0 

C(2F) 

33.1 


38.6 

39.1 

39.2 


Considering that errors of ~ 5 cm“^ in spin-orbit couplings are very small with 
respect to the error in the main non-relativistic energies we conclude that it is quite 
acceptable to compute spin-orbit couplings at a low level of theory. 


3. ScH molecule 


As discussed in the introduction and hinted by the results for the Sc atom presented 
in the previous section, from the point of view of high-resolution spectroscopy 
the accuracy expected for transition metal diatomics is much lower than the one 
generally achievable for molecules made up by main-group atoms. Also, because 
convergence seems to be rather irregular both with respect to the level of electron 
correlation and basis set size, one should not necessarily expect more expensive 
calculations to be much closer to experiment than simpler ones. 

For this reason, when possible, experimental data were used to adjust the ab 
initio potential energy curves, in particular, the experimental studies by Ram and 
Bernath 


19 


20(] where the u = 0 and sometimes v = 1 and v = 2 vibrational states 


of seven singlet terms (A^S+, A^A, R^II, and G'^II) were 

characterised. Of these the X, A and B terms dissociate to channel 1 of table [H C 
to channel 7 or perhaps 10, F to channel 6 while for terms D and G channels 3, 6, 
7, or 10 are all possible on the basis of symmetry considerations. 

Details on the refinement are given in section 01 in the rest of this section we 
discuss the ab initio calculations. 
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3.1. Potential energy curves 


Energy curves for the six molecular electronic terms correlating with the ground 
atomic states (dissociation channel 1 of tabled]) were computed using CASSCF and 
internally-contracted MRCI 37| in conjunction with the recent aug-cc-pwCVnZ 
basis sets (awcnz for short) [33, l3^- CASSCF orbitals (state-averaged over all the 
degenerate components of the six terms considered in this work) were used as a 
basis of the MRCI runs. 

A four-electron, ten-orbital complete active space comprising the scandium 
4s, 3d, 4p orbitals and the hydrogen Is orbital (5 active orbitals of ai symmetry, 2 
bi, 2 b 2 and 1 a 2 in the C 2 t; point group) was used in the calculations. The outer- 
core scandium 3s, 3p orbitals were left doubly occupied at the CASSCF stage but 
were correlated at the MRCI one. The inner-core Is, 2s, 2p orbitals were not corre¬ 
lated. As discussed by Balabanov and Peterson 3^, in multireference calculations 
the late transition metals Fe-Cu require an active space larger than the full-valence 
one, which should include a further set of diffuse d functions. However this is in 
not necessary for scandium. Inclusion of the Sc 4p orbitals is not thought to be 
indispensable for a correct description of bonding but was found to be necessary 
in practice to avoid convergence problems at the CASSCF stage. All curves were 
computed in the range 2.0 to 8.5 ao in steps of 0.05 ao and from 9.0 to 13.5 ao is 
steps of 0.5 ao, for a total of 141 points. 

Our best ab initio results are based on MRCI using the awc5z basis set; com¬ 
puting at this level the energies for all six terms for a single geometry takes about 
4 GB of RAM, 20 GB of disk space and 12 hours on a single core of an Intel Xeon 
E5-2670 CPU at 2.60 GHz. Potential energy curves include a relativistic correction 
computed as expectation value of the MVDl operator. The Davidson correction 
was not included in our final ab initio curves because, as already registed for the ex¬ 
citation energies of the scandium atom (see table [2|), it does not improve agreement 
with known experimental data; furthermore tests (not discussed here in detail) at 
the frozen-core / cc-pVDZ level showed that Davidson-corrected energies agreed 
worse with full Cl than uncorrected MRCI ones. 

The ah initio energy curves were slightly shifted in energy (i.e., their Eg were 
changed) so that they are exactly degenerate upon dissociation; in our MRCI cal¬ 
culation the exact degeneracy of the terms at r = -|-oo is broken mainly because 
the energies of (singlet or triplet) S"*" and A terms is computed simultaneously 
with a two-state calculation, while the (singlet or triplet) H terms were computed 
with one-state calculations; because of the internal contraction approximation used 
in Molpro in a multi-state calculation the variational flexibility of the MRCI wave 
function is increased and this results in a small downwards shift in energy. As a con¬ 
sequence at dissociation the two H terms are about 50 cm“^ higher in energy than 
the other terms; furthermore a small breaking of exact degeneracies is expected and 
normal also in uncontracted MRCI calculations because of the incomplete treat¬ 
ment of electron correlation. We therefore thought it was reasonable to shift all 
terms to restore the exact degeneracy. The shifts applied to MRCI/awc5z/MVDl 
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18000 



5 6 

bond length in bohrs 


Figure 1. Ab initio potential curves for ScH computed with MRCI and the awc5z basis set and the 
corresponding relativistic MVDl correction curves (see text for details). The 3s3p orbitals were correlated. 


curves for the 5^11, o ^A, 6^11 and terms are respectively 1.1, 49.7, 

1.5, 49.9 and 2.1 cm“^. The ground X term was taken as a reference and not 
shifted. 

Figured] presents our computed potentials. As it can be seen the ground X 
curves are distinct from the other curves: it has a much shorter equilibrium bond 
length and its relativistic correction curve is also very different from the others. This 
is a consequence of the different bonding character of this term discussed 

in section dl 

Equilibrium bond lengths rg, harmonic vibrational frequencies Wg and adiabatic 
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Table 4. Computed ab initio adiabatic electronic excitation energies (in cm ^), equilibrium bondlengths (in ao) 
and harmonic frequences (in cm~^) for selected ScH electronic terms. 

Term Anglada et al“ Hubert et al^ This work‘d 



Tg 

re 

Ule 

Tg 

re 

Ule 

Tg 

re 

Ule 

0 

3.41 

1621 

0 

3.35 

1611 

0 

3.34 

1587 


6600 

3.68 

1541 

5362 

3.58 

1439 

3914 

3.59 

1428 

B^U 

8400 

3.64 

1451 

— 

— 

— 

5856 

3.55 

1380 

a^A 

4600 

3.66 

1460 

3660 

3.55 

1450 

1868 

3.56 

1432 

b^n 

6200 

3.64 

1438 

— 

— 

— 

3544 

3.55 

1406 

c^E+ 

7900 

3.68 

1389 

— 

— 

— 

6122 

3.63 

1325 


“ Ref. 


13l |: values of the Tg’s are taken from the column labelled ‘B(lf)’ of table 8, 


re’s and_We’s from table 10. 
Ref, 


id ], using the data from the column labelled CCSDT12 (QC) of Table V 


and adding the relativistic corrections in the last column of Table VI. 

Using MRCI/awc5z/MVDl. The relative to state i was computed by 


uje'^ = \JV''{re’)/^ while the adiabatic excitation energy was computed as 
Vi{re^) — Vo{ri^^), where Vq is the potential for the ground state. 


„(0 


(0 


excitation energies Tg are reported in table S] and compared with previous theoret¬ 
ical calculations. 

Our computed equilibrium bond lengths are in very good agreement (within 
0.01 ao) with the recent theoretical values by Hubert et al 16 1. 

We compare our ab initio and empirically-refined results with energy levels re¬ 
constructed from the experimental study 20j in table [5l empirical refinement is dis¬ 
cussed in section 01 We prefer to compare directly with experimental energy levels 
because experimentally-deduced values for Tg, rg and Wg also include in an effective 
way spin-orbit and other coupling effects between different electronic terms. 

As discussed above, we expect our computed adiabatic excitation energies Tg to 
have errors of several hundreds or perhaps a few thousands of cm“^ and therefore 
they cannot be considered very accurate. On the other hand the equilibrium bond 
lengths and the shape of the potentials should be reasonably accurate, see also 
table O where the accuracy of the potential energy curves of the singlet states 
can be assessed by comparing to the vibrational energy separations within each 
state. The rotational intervals between the two lowest-J u = 0 energy levels are 
reproduced by our ab initio curves with errors of 0.00, 0.12 and 0.28 cm“^ for the 
X, A and B term respectively. Errors in vibrational transitions u = 0 to u = 0 are 
8.44 cm“^ for the X term and 29.27 cm“^ for the B term. 

For the ground X term only, we considered the error of the coupled-cluster based 
potential curves. The absolute errors for the u = 0 to u = 1 transition wavenumber 
using CCSD, CCSD(T), CCSDT and CCSDT(Q) are, respectively, 54.73, 24.45, 
6.75 and 1.15 cm“^ with the awc5z basis set, the DKH Hamiltonian and keeping all 
the coupling terms computed with MRCI or CASSCF; the CCSDT and CCSDT(Q) 
curves were obtained in the basis set formed by wc3z for Sc and 2z for hydrogen, and 
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Table 5. Selected energy levels for the lowest-energy singlet terms of ScH in cm~^. J is the total angular 
momentum (neglecting nuclear spin), v is the vibrational quantum num ber, p the parity; ‘Obs.’ are the v alue s 
derived using the spectroscopic constants reported by Ram and Bernath I 20 I 1 and the program PGOPHER Isdl : 
‘A’ and ‘R’ are the term values calculated with Duo using the ab initio (MRCI/awc5z/MVDl) and refined curves 
as described in the text. Calculations included spin-orbit and other couplings between electronic terms. 


J 

P 


energies 


obs - 

calc 



obs 

ab initio 

refined 

ab initio 

refined 




Xis+, 

V = 0 



0 

+ 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

- 

10.73 

10.73 

10.73 

0.00 

0.00 

10 

+ 

586.89 

586.88 

586.93 

-0.04 

-0.05 




Xis+, 

V = 1 



0 

+ 

1546.97 

1538.53 

1547.10 

8.44 

-0.12 

1 

- 

1557.45 

1549.00 

1557.54 

8.45 

-0.09 

10 

+ 

2120.06 

2111.48 

2118.31 

8.58 

1.76 




<r 

< 

u = 0 



2 

+ 

4213.36 

3830.27 

4213.35 

383.09 

0.01 

3 

+ 

4241.34 

3858.12 

4241.36 

383.21 

-0.03 

10 

+ 

4696.30 

4311.04 

4696.64 

385.25 

-0.34 

10 

- 

4696.30 

4311.01 

4696.61 

385.28 

-0.31 




B In, 

V = 0 



1 

+ 

5413.98 

5704.18 

5413.94 

-290.19 

0.04 

2 

+ 

5433.77 

5723.68 

5433.79 

-289.91 

-0.02 

10 

+ 

5943.02 

6226.10 

5939.67 

-283.07 

3.35 

10 

- 

5939.42 

6222.43 

5936.36 

-283.01 

3.06 




B In, 

V = 1 



1 

+ 

6776.75 

7037.67 

6776.70 

-260.92 

0.05 

2 

+ 

6795.96 

7056.69 

6796.35 

-260.73 

-0.39 

10 

+ 

7290.12 

7546.50 

7290.43 

-256.38 

-0.31 

10 

- 

7286.57 

7542.84 

7286.95 

-256.28 

-0.38 




B In, 

v = 2 



1 

+ 

8092.50 

8323.48 

8092.51 

-230.98 

-0.01 

2 

+ 

8111.15 

8342.00 

8111.07 

-230.85 

0.08 

10 

+ 

8590.66 

8818.77 

8589.10 

-228.10 

1.56 

10 

- 

8587.14 

8815.10 

8585.06 

-227.96 

2.08 


added as a correction to the CCSD(T) curve. These result indicate that our MRCI 
curves are, close to equilibrium, similar in quality to CCSDT and that quadruple 
excitations must be accounted for to obtain accuracies of the order of 1 cm“^. 


3.2. Dissociation energy 

The dissociation energy Dq of ScH is related to the potential well depth De by 

Do = De — Ezpe ( 1 ) 

where T^zpe = 787 cm“^ is the zero-point rotational-vibrational energy and the 
quoted value was computed using out MRCI/awc5z/MVDl PEC. The potential 
well depth Df, can be decomposed into three contributions: a main non-relativistic 
one, a spin-independent (scalar) relativistic contribution and a spin-dependent con¬ 
tribution due to spin-orbit: 


De = D 


NR 

e 


+ Df + D 


SO 

e 


( 2 ) 
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The spin-orbit contribution is due to the energy lowering of the scandium 
atom ^D 3/2 level with respect to the term and has a value = —3^/2 = 
— 101 cm“^, where ^ = 67.3 cm“^ is the atomic spin-orbit contant for the term 
(see table [I]). 

The (spin-orbit free) potential well depth computed with MRCI/awc5z is 
17459 cm“^; the Davidson correction gives a rather large shift of -1-1355 cm“^, 
leading for MRCI-|-Q/awc5z to a value D^^=18814 cm“^. 

With a view to ascertaining the quality of our ab initio curve close to dissociation 
we computed an accurate value for using high-order coupled cluster and the 
program MRCC 39]. 

Using the awc5z basis set and correlating the outer-core 3s3p orbitals gives for 
= 17694 cm“^ using CCSD and 18547 cm“^ using CCSD(T). The effect of full 
triples (T)—>-T was evaluated in a basis set formed by the wc3z for scandium and 
2z for hydrogen, giving a shift of -|-163 cm“^. The effect of quadruple excitations 
was evaluated in an even smaller basis set constucted complementing the 2z one for 
hydrogen and scandium with the core-correlation functions taken from the wc3z 
basis set and dropping the g functions; the computed shift T—>-Q is -|-48 cm“^; our 
best awc5z (frozen inner-core) coupled cluster value is = 18547 -|- 163 -|- 48 = 
18758 cm“^. The coupled cluster value therefore strongly supports the Davidson- 
corrected value for rather than the uncorrected MRCI one. 

We also considered the contribution to due to correlation of the inner-core 
2s2p orbitals; this effect gives a contribution of -|-62 cm“^ using CCSD/awc5z and 
-1-104 cm“^ using CCSD(T)/awc5z. 

The magnitude of basis set incompleteness was estimated by looking at the differ¬ 
ence between the awc5z values and the basis set extrapolated ones using the awc4z 
and awc5z basis sets; Basis set extrapolation of the awc5z values gives a very small 
contribution, namely —3 cm“^ for MRCI/awc5z and -|-6 cm“^ for CCSD(T)/awc5z. 

Finally, our best coupled-cluster-based value for the non-relativistic part of the 
potential well depth is = 18758 -|- 104 -|- 6 = 18868(50) cm“^, where the given 
uncertainty was assigned by halving of the sum of the quadruples correction, the 
difference between the CCSD and CCSD(T) inner-core corrections and the basis 
set extrapolation contribution. 

We now consider the scalar relativistic contribution D^. The MVDl/awc5z value 
is = -|-285 cm“^; using the DKH Hamiltonian (truncated to fourth order) and 
the awc4z-DK basis gives a contribution -|-311 cm“^ using MRCI, -1-291 cm“^ using 
MRCI+Q and +307 cm"! using CCSD(T). 

Taking the DKH value = 307 cm“^ (although there is no conclusive argument 
to favour it instead of the MVDl one) we arrive at a final value for the potential 
well depth Dg = 18868 + 307 — 101 = 19074(60) cm“^ and to a dissociation energy 
Dq = 19074 — 787 = 18287(60) cm“^, where the given uncertainty was increased 
to reflect the uncertainty on the relativistic correction. 


Kant and Moon 


40|| reported long ago an experimental value for the potential 
well depth Df, = 16613 + 700 cm“^ and Koseki et al 4]| gave a survey of calculated 
values of Dg which have a large spread of about 3000 cm“^ around the value quoted. 
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Our computed value is larger than the experimental one by about 2500 cm“^ and 
in disagreement with it by more than three times its uncertainty bar. 

Finally, we report some run-times for our coupled cluster results. The CCSDT 
run using the wc2z/2z basis set took for ScH 8.3 hours of CPU time on a 12-core 
Xeon X5660 at 2.80GHz machine (1.3 hours real time). The CCSDTQ run for ScH 
in the 2z-|-wC3z basis set took 10.6 days of CPU time (1.6 days real time) and 
5 GB of RAM on the same machine. A single GGSD(T)/awc5z run takes about 10 
minutes on a single core of the same machine. 


3.3. Dipole moment curves 


While potential energy curves can be refined semiempirically from (even limited) 
experimental data, one has normally to rely on computed, ab initio dipole mo¬ 
ment curves . With a view to computing accurate line intensities it is therefore 
important to produce dipole moment curves as accurate as possible. 


Le and Steimle 2l| reported for the ground term an experimental equi¬ 

librium dipole jjL = 1.74(0.15) D using optical Stark spectroscopy; for our work we 
choose the z axis such that a negative dipole corresponds to Sc'''H“ polarity, so we 
reverse their value to = —1.74(0.15) D. 

We were able to produce a very accurate value for the equilibrium dipole of the 
ground state term using coupled cluster and the energy-derivative (ED) method 
431] (A = ±10“^ au). As we are dealing with a closed-shell electronic state near 
equilibrium coupled cluster converges quickly with respect to the level of excitations 
included. High-order coupled cluster calculations used the program MRCG [3S|. 
Results are collected in table [6l our best theoretical value is pe = —1.72(2) D and 
is in full agreement with the experimental value of Le and Steimle 21 j. 

MRGI, on the other hand, has difficulties in reproducing the correct value for the 
dipole. Apart from the basis set size and whether or not core orbitals are correlated, 
we considered three further factors affecting MRGI dipoles, namely: i) whether 
dipoles are computed by expectation value (XP) or energy derivative (ED) [43l |: 
a) whether at the CASSCF step orbitals are obtained by state averaging (state- 
averaged orbitals, SAG) or are specifically optimised for for the electronic 

term (state-specific orbitals, SSO); and, iii)^ the effect of using Davidson-corrected 
energies instead of MRGI ones. We compared dipoles obtained with MRGI in the 
awc3z basis set with the very accurate (non-relativistic, 3s3p correlated) value ob¬ 
tained with coupled cluster, which gives in this basis set the value p, = 1.66 D 
(see table O GGSD(T)/awc3z value -|- lines B,G,D). Results are collected in ta¬ 
ble [71 Figure [2] shows the dipole moment curves for the X ground term computed 
with GCSD(T), MRGI/XP, MRCI/ED and MRGI-hQ/ED; note that the GGSD(T) 
curve become unphysical at large bondlengths. 

As one can see in table [ 7 ] both GASSGF and MRGI/XP equilibrium dipoles are 
too small in magnitude by about 0.3-0.5 D, a considerable amount; using state- 
specific orbitals reduces the error to about 0.3 D, indicating that the GASSGF and 
MRGI wave functions are quite far from the exact, full GI one (which is independent 
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Table 6. Equilibrium dipole of the ground state X term using coupled cluster theory. Dipoles were computed 
at r = 3.35 ao using the energy-derivative method and A = ±10~^ au. A part from the line labelled ‘D’, all 
calculations correlated the outer-core 3s3p orbitals but the kept the inner core Is2s2p uncorrelated. Dipoles are 
in debyes. The experimental value is -1.74(0.15) D I 21 I 1 . 


label 

method 

basis set 

value 


RHF 

awc3z 

-1.436 


CCSD 

awc3z 

-1.668 


CCSD(T) 

awc3z 

-1.640 


CCSD(T) 

awc4z 

-1.686 


CCSD(T) 

awc5z 

-1.702 

A 

CCSD(T) 

awe [345] z“ 

-1.719 


higher 

order correlation 

B 

(T) ^ T 

wc3z/2z^ 

-0.002 

C 

T ^ T(Q) 

wc3z/2z^ 

-0.009 

D 

T(Q) ^ Q 

2z-|-wC(3z)'^ 

-0.005 


inner- 

core correlation*^ 


E 

CCSD(T) 

awc3z 

-0.011 


relativistic'^ 


F 

CCSD(T) 

awc4z 

-h0.078 


vibrational averaging-^ 


G 



-0.054 

A + B +C +D +E+F+G 

best ab initio^ 


-1.72(2) 


“ Basis-set extrapolated value using a fin = fJ-e + A/n^ formula. 

^ Correction due to full triples and perturbative quadruple excitations using the 
cc-pVDZ basis set for hydrogen and cc-pwCVTZ for scandium. 

Correction due to full quadruple excitations using the cc-pVDZ basis set for hy¬ 
drogen and the cc-pVDZ complemented with the core-valence correlation functions 
from the cc-pwCVTZ (g functions excluded) for scandium. 

Correction due to correlation of the 2s2p orbitals. The innermost Is orbital was 
not correlated. 

® Correction due to scalar-relativistic effects computed as difference of 
CCSD(T)/DKH4/awc5z-dk and CCSD(T)/awc5z dipoles. 

^ Vibrational averaging computed as (0|/x(r)|0) — fi{r = 3.35ao), where |0) is the 
J = 0,u = 0 vibrational ground state (obtained using the ah initio MRCI/awc5z 
PEC) and fi{r) is the MRCI/energy-derivative/awc5z dipole moment curve. 

3 The estimated uncertainty in the teoretical dipole is mostly due to residual ba¬ 
sis set incompleteness error and incomplete treatment of higher-order correlation 
effects. 
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Table 7. Values and errors in computed MRCI dipoles at r = 3.35 ao for the ground state term. The 

acronym SAO and SSO stand for state-averaged orbitals and state-specific orbitals, respectively. The last two 
columns report differences with the accurate value ^ — 1.66 D obtained with coupled cluster (see text). Dipoles 
are in debyes. 


method^ 

values 

values 

- exact 


SAO 

SSO 

SAO 

SSO 

CASSCE/XP 

-0.99 

-1.32 

0.67 

0.34 

MRCI/XP 

-1.19 

-1.35 

0.47 

0.31 

MRCI/ED 

-1.56 

-1.66 

0.08 

0.00 

MRCI-hQ/ED 

-1.68 

-1.71 

-0.02 

-0.05 


“ XP or ED specifies whether dipole were computed as expectation value or energy 
derivatives (field strength A = ±10““^ an). 



Figure 2. Ab initio dipole curves for the ground X term computed in the awc5z basis set and various 
methods (see text). 


on the choice of the orbitals). Computing the dipole by energy derivative greatly 
reduces the error in the MRCI dipoles, bringing them in much closer agreement 
with the coupled cluster value. Using Davidson-corrected energies (MRCI-I-Q) in¬ 
troduces a shift of about 0.12 D and brings the dipole curve near equilibrium even 
closer to the coupled cluster one (see also hg. ED. The Davidson-corrected dipole 
curve (‘relaxed’ reference energies were used) has a small jump discontinuity (0.02 D 
of magnitude) at r = 4.6 ap; on the other hand the MRCI/ED curve is perfectly 
smooth for all bond lengths. In conclusion the major factor affecting MRCI dipoles 
is whether the XP or ED technique is used, with ED producing considerably better 
dipoles. Using SSO instead of SSA helps somewhat but is of secondary importance. 

It is also worth noting that, although the absolute value of MRCI/XP dipoles 
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Figure 3. Diagonal (in bold) and off-diagonal dipole moment curves for ScH computed with with MRCI 
and the awc5z basis set. 


is considerably off, this quantity only affects line intensities for pure rotational 
transitions. Intensities of vibrational transition within an electronic term depend 
on the shape of the dipole function, and may be given quite accurately even by 
MRCI/XP. 

Finally, we decided to use awc5z/MRCI/ED for all diagonal dipole curves in 
virtue of their smoothness, although using MRCI+Q dipoles may lead to a slight 
improvement. Off-diagonal dipoles were computed as expectation values of the 
awc5z/MRCI wave functions; although it is possible to compute off-diagonal dipoles 
using an energy-derivative technique 4^, this route was not pursued at this time. 
Figure E] shows the diagonal and off-diagonal dipole moment curves for the various 
electronic. 


3.4. Spin-orbit and other coupling curves 


We computed spin-orbit couplings and couplings of the angular momentum oper¬ 
ators Lx and Ly using the CASSCF or MRCI wave functions. 

Figure 0] shows matrix elements of the Lx and Ly operator, obtained at the 
CASSCF/awc3z level; these couplin gs e nter in the L-uncoupling and spin-electronic 
terms of the rotational Hamiltonian 4^ and are responsible for A-doubling. Finally, 


fig. [5] reports the 10 symmetry-independent spin-orbit coupling curves obtained 
at the CASSCF/awc3z level. Care was taken to ensure that the coupling curves 
and dipoles are both smooth and phase corrected, something that is by no means 
standard in the literature 


46 |. 
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Figure 4. Matrix elements of the Lx and Ly operators for ScH computed with CASSCF and the awc3z 
basis set. Specifically, the curve labelled ‘1’ is IS.^2_y2\Lx\^'^y)/i\ curve ‘2’ is (}-/\^2_y2\Lx\^'^y)/(—i)\ 
curve ‘3’ is {^E'*"|Lypna;}/(\/32) ; curve ‘4’ is {^S'*'|Ly|^na;}/(\/3z). The phases of the electronic wave 
functions are chosen such as (^I\.x\Lz\^L\.y) = 2 , (^TixlLzl^Tiy) = i, A^2_y2\Lz\^Axy) = —2i and 

A^2_y2|T2pAa3y} = 2i. 



Figure 5. Spin-orbit coupling matrix elements for ScH computed with CASSCF and the awc3z ba¬ 
sis set. Specifically: the curve labelled ‘F is {^A^ 2 _y 2 ,S = llH’goPAccy, S = l}/(—i); curve ‘2’ is 
Aa;y|^So|^Aa;x-yy, s = 0}/(— 2 ); curve ‘3’ is {^S+,S = 0|H'so|^ny,S = l}/(— 2 ); curve ‘4’ is 
{iny|J^soPS+,S = l}/ 2 ; curve ‘5’ is (iS+lJ^soPn^, E = l>/ 2 ; curve ‘6’ is {3n^,S = l|Hsopny,E = 
l}/(— 2 ); curve ‘7’ is {^Hxl-HsoPriy, S = 0}/2; curve ‘8’ is {^nx,E = OlHgoP Aa:^, E = 1)/^; curve ‘9’ is 
Axy I-Hgopna;, E = 1}/(— 2 ); curve ‘10’ is {^Hy I^Hgol^A^a;—yyi ^ — 1}/(~0- The phases of the electronic 
wave functions are chosen such as (^Ha^lLzI^Hy) = i, (^Tix\Lz\^Tly) = 2, A^2_y2|L2|^Aa;y} = —22 and 

A^2_y2 |Z/z|^ Axy) = 22. 
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4. Line list 


The potential energy, dipole and coupling curves were then used with the in-house 
program Duo to produce a line list for ^^ScH. Duo solves in an essentially exact 
way the rotational-vibrational-electronic problem for multiple interacting energy 
curves for diatomic molecules and is described in detail elsewhere 28l. l47l. 1481] . The 


line list can be obtained from www.exomol.com, while all the curves used to produce 
it are made available as supplementary material. In all nuclear-motion calculations 
we used the atomic masses mn = 1.0078250321 Da and msc = 44.9559100 Da, 

= 0.985726930 Da = 


-|- m 


Sc 


which give for ScH a reduced mass // = (rujj 
1796.87027 me. 

The ab initio potential energy curves of the singlet terms were adjusted by fitting 
to the energy term values (J < 12) derived from the experimental spectroscopic 
constants reported by Ram and Bernarth 20|] which cover vibrational excitations 
with u = 0,1 (X), V = 0 (^) and u = 0,1, 2 (B) only. This was not possible for 
the triplet states (see table SD- We also decided not to use in the adjustments the 
very recent experimental data on the electronic state 22] since it lead to an 
equilibrium bond length Cg = 3.94 ao which differs too substantially from theory 
to be safely trusted, see table [H 

In the refined curves the dissociation energy was fixed to the Dg value by Koseki 
et al [^. The triplet curves were also scaled, to dissociate to the same value of 
Dg. All refined curves are given as supplementary material to the paper together 
with the ah initio curves. The triplet electronic states appear to be in the strong 
resonance with the rovibronic states from B and prevented an accurate fit to 
the R-state energies, especially for J > 12. It should be noted however that the 
spectroscopic constants from Ref. 20(| were derived using in the absence of the of 
interaction with the triplet states, and thus can also contain artifacts. 

We then used the program Duo [4^ to solve the coupled Schrodinger equations 
to compute the rovibronic energies of ScH up to the dissociation. In particular, we 
obtained for ^^ScH a zero-point-energy of 799.6 cm“^. The highest value the total 
angular momentum J can assume for bound states is found to be J = 59. 

The corresponding rovibronic eigenfunctions were combined with the ah initio 
dipole moment curves to produce Einstein A coefficients for all transitions with 
line positions up to Dq. The Einstein A coefficients together with the rovibrational 
energies supplemented by the total degeneracies and quantum numbers make up 
the line list. 

Duo calculations consist of two steps: in the first step we used a grid of 501 points 
to solve six separate vibrational Schrodinger equations for each electronic state, 
using as potential the ah initio potential curves shown in Eig. [T]or the empirically 
adjusted curves described above. We then selected 40 lowest-energy eigenfunctions 
from each set; the union of these 40 x 6 = 240 functions constitutes our vibrational 
basis set |state,u), where v is the vibrational quantum number and ‘state’ is the 
label identifying the electronic state. In the second step of the calculation we build 
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a basis set of Hund’s case a functions of the type 


|state, A, 5, S, u) = |state, A, 5, S, )| J, U, M)|state, u), 


( 3 ) 


where |state, A, S, S) is the electronic function, |J, n,M) is the rotational func¬ 
tion and |state, u) is one of the vibrational functions; A, S, and are the z axis 
projections of the electronic, spin and total angular momenta, respectively, and 
n = A -|- S; M is the projection of the total angular momentum along the labora¬ 
tory axis Z. The full, coupled problem is then solved by exact diagonalization in 
the chosen basis set. 

In order to guarantee that all phases of the ab initio couplings as well as transition 
dipole moments are consistent, we used the matrix elements of the operator 
between the corresponding degenerate 11 and A components as provided by Molpro. 
These matrix elements were then used to transform the matrix elements of all 
coupling to the representation where Lz is diagonal, which is used in Duo according 
with eq. ([3]). The phases are chosen such that all matrix elements are positive. 

Our ^^ScH line list contains 1 152 827 transitions and is given in the ExoMol 
format 49|] consisting of two files, an energy file and a transition file. This is 


based on a method originally developed for the BT2 line list Extracts for the 
line lists are given in Tables [8] and [9l Using all energies of ^^ScH we computed 
its partition function for temperatures up to 5000 K. The line lists and partition 
function together with auxiliary data including the potential energy, spin-orbit, 
electronic angular momentum, and dipole moment curves, as well as the absorption 
spectrum given in cross section format 5l|, can all be obtained from the ExoMol 
website at www.exomol.com. 

As an example, in fig. [6] we show absorption (T = 298 K) intensities of ScH 
generated using a stick-diagram. Figure [3 illustrates in detail the band structure 
of the absorption spectrum of ScH at T = 1 500 K generated using a Gaussian line 
profile with the half-width-at-half-maximum (HWHM) of 5 cm“^. As one can see, 
the strongest features belong to the X-X, B-X, a-a, and c-b electronic bands. 
The triplet bands electronic bands a-a and c-b should be strong enough to be 
potentially observable in the lab. Due to the spin-orbit couplings between different 
components forbidden bands also appear to contribute, with AX and cX being the 
strongest. These bands are significantly weaker than the dipole-allowed ones, but 
could still represent an important source of ScH opacity. 


5. Conclusions 

A hot line list containing pure rotational, ro-vibrational and vibronic transitions 
for ScH was generated using new ab initio potential energy, dipole moment, spin- 
orbit, and electronic angular moment curves obtained at a high level of theory. The 
work was performed with a view to astrophysical applications. The analysis of the 
importance of different absorption bands for the opacities of ScH is presented. 
The complexity of the electronic structure problem when transition metals are 
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Table 8. Sample extract from the energy file for "^^ScH. The whole file contains 8 451 entries. 


i 

E 

9 

J 

+/- 

e// 

State 

V 

|A| 

|S| 

IQI 

1 

0.000000 

16 

0 

+ 

f 

XlSigma+ 

0 

0 

0 

0 

2 

1547.095548 

16 

0 

+ 

f 

XlSigma+ 

1 

0 

0 

0 

3 

3019.322250 

16 

0 

+ 

f 

XlSigma+ 

2 

0 

0 

0 

4 

3352.480112 

16 

0 

+ 

f 

b3Pi 

0 

1 

1 

0 

5 

4430.504406 

16 

0 

+ 

f 

XlSigma+ 

3 

0 

0 

0 

6 

4707.209870 

16 

0 

+ 

f 

b3Pi 

1 

1 

1 

0 

7 

5789.270187 

16 

0 

+ 

f 

XlSigma+ 

4 

0 

0 

0 

8 

6015.690883 

16 

0 

+ 

f 

b3Pi 

2 

1 

1 

0 

9 

7100.259438 

16 

0 

+ 

f 

XlSigma+ 

5 

0 

0 

0 

10 

7277.938899 

16 

0 

+ 

f 

b3Pi 

3 

1 

1 

0 

11 

8363.776161 

16 

0 

+ 

f 

XlSigmat- 

6 

0 

0 

0 

12 

8494.401968 

16 

0 

+ 

f 

b3Pi 

4 

1 

1 

0 

13 

9574.352490 

16 

0 

+ 

f 

XlSigina+ 

7 

0 

0 

0 

14 

9667.692015 

16 

0 

+ 

f 

b3Pi 

5 

1 

1 

0 

15 

10718.414759 

16 

0 

+ 

f 

b3Pi 

6 

1 

1 

0 

16 

10804.884934 

16 

0 

+ 

f 

XlSigina+ 

8 

0 

0 

0 

17 

11788.395483 

16 

0 

+ 

f 

b3Pi 

7 

1 

1 

0 

18 

11902.873682 

16 

0 

+ 

f 

XlSigma+ 

9 

0 

0 

0 

19 

12788.251034 

16 

0 

+ 

f 

b3Pi 

8 

1 

1 

0 

20 

12941.194344 

16 

0 

+ 

f 

XlSigma+ 

10 

0 

0 

0 

21 

13714.254989 

16 

0 

+ 

f 

b3Pi 

9 

1 

1 

0 

22 

13898.140315 

16 

0 

+ 

f 

XlSigma+ 

11 

0 

0 

0 

23 

14551.459500 

16 

0 

+ 

f 

b3Pi 

10 

1 

1 

0 

24 

14746.814942 

16 

0 

+ 

f 

XlSigma+ 

12 

0 

0 

0 


i: State counting number. 

E: State energy in cm“^. 
g: State degeneracy. 

+/—: Actual state parity. 

e/f: Rotationless parity. 

v: State vibrational quantum number. 

|A|: Absolute value of A (projection of the electronic angular momenum). 

|S|: Absolute value of E (projection of the electronic spin). 

|Q|: Absolute value of = A + E (projection of the total angular momentum). 


Table 9. Sample extract from the transition file for “^^ScH. The whole file contains 1 152 826 entries. 


/ 

i 

Aif 

1351 

1231 

3.3006E-07 

3574 

3468 

3.8843E-07 

5782 

5693 

5.0269E-06 

4942 

5037 

9.9272E-06 

7688 

7734 

4.9607E-03 

2070 

1952 

3.8782E-01 

2919 

2580 

1.0196E+00 

2804 

2692 

1.0196E+00 

1362 

1713 

8.5042E-08 

5638 

5727 

6.8009E-04 

3137 

3242 

4.2340E-06 

3672 

3568 

1.9332E-04 

7406 

7350 

4.6729E-06 

5984 

5901 

1.2467E-05 

3396 

3505 

9.4844E-05 

993 

1110 

1.0904E-06 

3324 

2999 

1.9023E-05 

2398 

2282 

1.0987E-04 


/: Upper (final) state counting number. 
i: Lower (initial) state counting number. 
Ajf: Einstein A coefficient in s“^. 
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Figure 6. Absorption cross-sections of ScH at T=298 K. 
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Figure 7. Overview of the absorption line intensities of ScH at T=1500 K. 


involved means that the accuracy of these calculations is much worse than what is 
normally achievable for small molecules containing light main-group elements. To 
help mitigate this problem, it is desirable to utilise experimental data to improve 
the accuracy of the results. We have done this for the singlet states. 

However, for transitions involving triplet states there are no measured data for 
us to use: for example, the recent triplet-triplet measurements by Mukund et al 
22l | in the 17 940 cm“^ region lie well above our calculated transitions and have 


not reported any absolute energies of the lower a^A electronic state. Besides, their 
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only parameter with which we can directly compare, appears to be too low 
to be fully trusted. Thus for the triplet states we used the Tg valued taken from 
our ah initio calculations, which may be in error by up to a few thousands cm“^, 
implying that the entire bands may be erroneously shifted by similar amounts. 
Also, structure due to perturbations in both the singlet and triplet manifolds are 
dependent on the potential energy curve separations and therefore may not be 
reproduced accurately. However we expect that our line list to be quantitatively 
accurate enough for the singlet state energies and to provide a detailed spectral 
structure of individual bands to be useful. 
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